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We review the most recent results, derived within the combined framework 
of the pinch technique and the background field method, describing certain 
QCD nonperturbative properties. 

Since its introduction, the combined framework of the pinch technique 
(PT) [U [21 [3j HI [5] and the background field method [6] , known in the lit- 
erature as the PT-BFM scheme [3 El [9] , has provided a sound theoretical 
basis for addressing the nonperturbative study of the QCD Green's func- 
tions of both the gluon as well as the ghost sector of QCD, respecting at 
the same time the fundamental symmetries of the theory. Results derived 
within this framework include, but are not limited to, the first evidence 
of the existence (in the Landau gauge) of massive solutions in the prop- 
erly truncated QCD Schwinger-Dyson equations (SDEs) as found in lattice 
simulations [10] - and which can be interpreted in terms of a nonpertur- 
bative mass [TQ1 [HI 02] which tames the infrared (IR) divergences of the 
Green's functions of the theory-, the study of the Kugo-Ojima function [T3] 
and the identification of the role of the ghost for achieving a chiral symme- 
try breaking pattern that provides for dynamically generated quark masses 
compatible with phenomenology [14J. 

In this talk I will present the latest results derived within the PT-BFM 
framework (in the Landau gauge) and discuss in particular: 

• The use of the SDEs to compute the nonperturbative modifications 
caused to the IR finite gluon propagator by the inclusion of a small 
number of quark families [15J; 
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• The general derivation of the full non-perturbative equation that gov- 
erns the momentum evolution of the dynamically generated gluon 
mass [IB] . 



As described in [15] the PT-BFM allows to develop an approximate 
method for "unquenching" the (IR finite) gluon propagator, computing non- 
perturbatively the effects induced by a small number of light quark families. 
The procedure consists of two basic steps: 

• Computing the fully-dressed quark-loop diagram, using as input the 
nonperturbative quark propagators obtained from the solution of the 
gap equation, together with an Ansatz for the fully-dressed quark- 
gluon vertex that preserves gauge-invariance |14j ; 

• Adding this result to the quenched gluon propagator obtained in large- 
volume lattice simulations. 

The key assumption of the method sketched above, is therefore that the 
effects of a small number of quark families to the gluon propagator may be 
considered as a "perturbation" to the quenched case, of which the quark- 
loop diagram constitutes the leading correction term, with the subleading 
terms stemming from the (originally) pure Yang-Mills diagrams which now 
get modified from the quark loops nested inside them. Thus, within the 
approximations we will employ these latter corrections are neglected, so 
that one can identify (even when dynamical quarks are present) with the 
quenched lattice propagator all SDE graphs except the quark loop diagram. 

The expression for the PT-BFM scalar cofactor A Q (g 2 ) of the unquenched 
propagator (the subindex "Q" standing for "quarks"), defined as 
&Q U (q) = P M y(g)A Q (g 2 ) with = - q^q u /q 2 the dimensionless trans- 
verse projector, can be then written as [15] 



In what follows we will describe all the different terms appearing in the 
right-hand side of the formula above. 



1. Unquenching the gluon propagator 



A Q (g 2 ) 



A(9 2 ) 



(1) 



1 + 



{iX(q2)[l + G(q 2 )}- 2 -A 2 }a(? 2 ) 



A(g 2 ) is the quenched propagator which, as already pointed out, will 
be identified with the one obtained from the large volume lattice sim- 
ulations. 
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• G(q 2 ) is a special Green's function particular to the PT-BFM which 
achieves the conversion from the PT-BFM to the conventional gluon 
propagator [3 [9] ; in the Landau gauge it is known to coincide with 
the Kugo-Ojima function [131 [T7]. 

• X(q 2 ) is the PT-BFM scalar cofactor resulting from the calculation of 
the quark loop diagram; defining X^^q) = P flu (q)X(q 2 ), one has 



where S is the full fermion propagator (with S 1 (p) = —iA(p)[p l — 
M(p)] and A4 the dynamical quark mass), and T v is the full PT- 
BFM quark-gluon vertex. To evaluate expression ([2]) one proceeds as 
follows: (?) the nonperturbative behavior of the functions A and M 
appearing in the definition of the full quark propagator are obtained by 
solving numerically the quark gap equation as done in |14j ; (ii) for the 
full PT-BFM quark-gluon vertex T one uses a suitable nonperturbative 
Ansatz, satisfying the gauge symmetry of the theory -such as the Ball- 
Chiu vertex [18J or the Curtis-Pennington vertex [19]- improved with 
the inclusion of the (numerically crucial) dependence on the ghost 
dressing function and the quark-ghost scattering amplitude [14 . To 
be sure, other forms of the quark-gluon vertex exists, such as those 
reported in [2i)\ 121] . and it would be interesting to check what effects 
they might have on our predictions. 

• Finally, A 2 = Ag X (0) - A _:L (0) = m|(0) - m 2 (0) denotes the gluon 

mass difference at q 2 = (notice that since X(0) = 0, the quark 
contribution to this quantity is only indirect,! e., through the modifi- 
cation it will induce on the various ingredients appearing in the mass 
equation - see next section) . A solid first-principle determination of A 2 
has not been attempted, mainly due to the fact that the derivation of 
the complete mass equation has been only very recently achieved |16] 
(see the next section again); in the analysis presented here we will 
restrict ourselves to extracting an approximate range for A 2 , by em- 
ploying a suitable extrapolation of the (unquenched) curves obtained 
from intermediate momenta towards the deep IR. 

The main results of our study may be summarized as follows (see Fig. [T|) . 
The basic effect of the quark loop(s) (one or two families with a constituent 
mass of the order of 300 MeV) is to suppress considerably the gluon prop- 
agator in the IR and intermediate momenta regions, while the ultraviolet 
tails increase, exactly as expected from the standard renormalization group 
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Fig. 1: (Top panels) Solution of the quark gap equation: A~ 1 (k 2 ) (left) and 
dynamical quark mass M.(k 2 ) (right) renormalized at \i = 4.3 GeV; dotted 
black curves correspond to the improved Ball-Chiu vertex, while dashed blue 
curves to the improved CP vertex. (Bottom panels) Comparison between the 
quenched and the unquenched gluon propagator (left) and dressing function 
(right). The yellow striped band in the left plot shows the possible values 
that Aq(0) can assume at zero momentum depending on the extrapolation 
point used; in the case of the dressing function (which is basically insensible 
to the IR saturation point) we used a curve with an extrapolation point at 
q 2 = 0.05 GeV 2 . The quenched lattice results of [22] are also displayed for 
comparison. 



analysis. In addition, the inclusion of light quarks makes the gluon propa- 
gator saturate at a lower point, which can be translated into having a larger 
gluon mass. As far as the gluon dressing function Z(q 2 ) = q 2 A(q 2 ) is con- 
cerned, one observes a suppression of the intermediate momentum region 
peak. A comparison with the recent full QCD lattice simulations of [23] is 
currently underway; however a comparison with some of the available lat- 
tice data [23] (Fig. [2]) shows an excellent qualitative agreement as well as 
a rather favorable quantitative agreement (with discrepancies at the 20% 
level maximum). 



Binosi-eQCD2012 printed on August 30, 2012 



■5 







Gluon propagator 






* Lattice n =2 and V= 1 6 3 X 36 






* Lattice n =2 and V= 12 3 X 24 






a Lattice n=0 and V= 64 4 






* Lattice n =0 and V= 72 4 






* Lattice n=0 and V= 80 4 






n=2and 1^=710 MeV 






n=2 and M^S MeV 






fitforn=0 















q [GeV] 





Gluon dressing function 




+ Lattice n=2and V=16 ! x32 




h Lattice n =0 and V=64 4 




s> Lattice n=0andV=72'' 


ly-A 


4- Lattice n =0 and V=80 4 




n=2 and M m =7A0 MeV 


" I • I%\ 


n=2and 1^=15 MeV 




fitforn=0 







q[GeV| 



Fig. 2: The unquenched gluon propagator (left) and dressing function (right) 
obtained in [24J (dark gray stars), together with the SDE results for two light 
quarks with M u j d = 15 MeV (dashed blue curve ) and M u / d = 710 MeV 
(dotted black curve). Quenched data |22| are again shown for comparison. 



2. The complete gluon mass equation 

Massive solutions of the gluon propagator SDE can be parametrized 
as (Euclidean space) A~ 1 (q 2 ) = q 2 J(q 2 ) + m?(q 2 ); therefore one faces the 
fundamental question of how to disentangle from the SDE the part that 
determines the evolution of the mass m 2 (q 2 ) from the part that controls 
the evolution of the "kinetic" term J(q 2 ). This is to be contrasted to what 
happens in the analogous studies of chiral symmetry breaking, where one 
derives a system of two coupled equations, one determining the "wave func- 
tion" ("kinetic part") of the quark self-energy, and one determining the 
dynamical (constituent) quark mass [141 125j . Of course, in the case of the 
quark self-energy the above separation of both sides of the corresponding 
SDE (quark gap equation) is realized in a direct way, due to the distinct 
Dirac properties of the two quantities appearing in it, while in the case of the 
gluon propagator no such straightforward separation is possible. However, 
an unambiguous way for implementing this separation, which exploited to 
the fullest the characteristic structure of a certain type of vertices that are 
inextricably connected with the process of gluon mass generation and nat- 
urally appears in the PT-BFM framework, was recently presented in [16] . 

Specifically, a crucial condition for obtaining out of the SDEs an IR- 
finite gluon propagator without interfering with the gauge invariance of the 
theory, is the existence of a set of special vertices that are purely longitu- 
dinal and contain massless poles, and must be added to the usual (fully- 
dressed) vertices of the theory. The role of these vertices is two-fold. On 
the one hand, thanks to the massless poles they contain, they make possible 
the emergence of a IR finite solution out of the SDE governing the gluon 
propagator; thus, one invokes essentially a non-Abelian realization of the 
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well-known Schwinger mechanism [26} [27] . On the other hand, these same 
poles act like composite Nambu-Goldstone excitations, preserving the form 
of the STIs of the theory in the presence of a gluon mass. 

It turns out that the very nature of these vertices furnishes a solid guid- 
ing principle for implementing the aforementioned separation between mass 
and kinetic terms. In particular, their longitudinal structure, coupled to the 
fact that one works in the Landau gauge, completely determines the longi- 
tudinal component of the mass equation; this is tantamount to knowing the 
full mass equation, given that the answer is bound to be transverse. 

Due to the complexity of the derivation of the equation, we will not 
discuss it here, but rather sketch its final form as well as its main ingredients, 
together with the numerical solutions it gives rise to. 

Schematically the equation reads 



where /Ci is the contribution coming form the one-loop dressed diagrams 
(namely the graphs appearing in the PT-BFM gluon propagator SDE con- 
taining trilinear vertices only), whereas /C2 is the contribution of two-loop 
dressed diagrams (that is, the graphs containing quadrilinear vertices). As 
indicated in Eq. ([3]), while IC± contains only the gluon propagator, in /C2 a 
new form factor Y appears, which involves the three gluon vertex and reads 



The lowest order perturbative calculation of Y (obtained by substituting 
tree-level values for all quantities appearing in the expression above) yields 
(after renormalization) Y ~ log/c 2 //U 2 ; this value multiplied by a constant 
C (basically modelling, in a rather heuristic way, further corrections that 
may be added to the "skeleton" provided by the lowest order result) is the 
one used in [16] for studying numerically the solutions of Eq. Q. Of course, 
the value of C corresponding to the lowest order expression is fixed to the 
actual value C = 3ttCaO( s ; however, it is convenient to treat C as a free 
parameter, thus disentangling it from the value of a s , and studying what 
happens to the solution spectrum of Eq. ([3]) when the two parameters are 
varied independently. 

As shown in the left panel of Fig. [3] (where C is now measured in units 
of 3itCa), there is a continuous curve formed by the pairs (C, a s ), for which 
one finds physical solutions. Indeed, for small values of C one has that no 
solution exists; this absence of solutions persists (for the quenched case) 
until the critical value C ~ 0.56 is reached, after which one finds exactly 
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Fig. 3: The curve described by the set of the pairs (C,a s ) for which one 
finds physical solutions to the full mass equation ([3]) (left), and the typical 
monotonically decreasing solution of the mass equation @ (right). 



one monotonically decreasing solution. However, for values up to C ~ 0.8 
the coupling needed to get the corresponding running mass is of 0(1), while 
for the quenched case the expected coupling from the 4-loop (momentum 
subtraction) calculation is a s = 0.22 at = 4.3 GeV |28j. This latter 
value is obtained for C ~ 1.55 - 1.65, whereas for C ~ 0.88 one finds the 
solution to Eq. ([3]) for the lowest order perturbative value of the coefficient. 
In general one observes, as expected, that as C is increased, a s decreases, 
e.g., for C = 1.1, 1.3 and 1.5 one obtains solutions corresponding to the 
strong coupling values a s 0.51, 0.35 and 0.26, respectively. 

In the right panel of Fig. [3] we plot the solutions for the most represen- 
tative C values, i.e., C = 0.88 and C = 1.6 (corresponding to, as already 
said, a s ~ 0.88 and 0.22 respectively), normalized in such a way that the 
mass at zero coincides with the IR saturating value found in lattice (Landau 
gauge) quenched simulations [22J, or m 2 (0) = A _1 (0) 0.141 GeV 2 . As 
can be readily appreciated, the masses obtained display the basic qualitative 
features expected on general field-theoretic considerations and employed in 
numerous phenomenological studies; in particular, they are monotonically 
decreasing functions of the momentum, and vanish rather rapidly in the 
ultraviolet [DIMIED]- It would seem, therefore, that the PT-BFM all-order 
analysis described here puts the entire concept of the gluon mass, and a 
variety of fundamental properties ascribed to it, on a solid first-principle 
basis. 
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